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Abstract 

We prove that the focusing cubic wave equation in three spatial dimensions 
has a countable family of self-similar solutions which are smooth inside 
the past light cone of the singularity. These solutions are labeled by an 
integer index n which counts the number of oscillations of the solution. 
The linearized operator around the n-th solution is shown to have n + 1 
negative eigenvalues (one of which corresponds to the gauge mode) which 
implies that all n > solutions are unstable. It is also shown that all n > 
solutions have a singularity outside the past light cone which casts doubt 
on whether these solutions may participate in the Cauchy evolution, even 
for non-generic initial data. 



1 Introduction 



This paper is a continuation of our studies of semihnear wave equations in three 
spatial dimensions with a focusing power nonlinearity 

dttu — Au — = , p = odd integer > 3 . (1) 

In [1] we showed that for each odd integer p > 7 equation ([T]) has a countable 
sequence of regular self-similar solutions while for p = 5 there is no nontriv- 
ial regular self-similar solution. This result has important consequences for the 
character of the threshold of blowup for Eq. ([1]) [2] . 

Here we consider the subcritical power p = 3 and, as before, restrict our 
attention to spherically symmetric solutions, so u = u{t, r) and Eq. ([T]) reduces 
to 

2 

dttU — drrU drU — = . (2) 

r 

This equation has the scaling symmetry (for each positive constant A) 

uit, r) ^ uxit, r) = X-\{t/X, r/X) , (3) 

so it is natural to ask if there are solutions which are invariant (modulo time 
translation) under this scaling. Such solutions are called self-similar. It follows 
from ([3]) and the symmetry under time translation that self-similar solutions must 
have the form 

«(t,r) = (r-t)-if/(p), (4) 

where T is a positive constant (usually referred to as the blowup time) and 
p = r/{T — t) is the so-called similarity variable ranging from zero to infinity. By 
definition, the self-similar solutions are singular at the point (T, 0). 

Substituting the ansatz (jlj) into Eq. ([2]) we get the ordinary differential equa- 
tion 

The obvious constant solution of this equation is Uo{p) = and the question 
is if there are other nontrivial solutions U{p) which are smooth in the interval 
< p < 1 (i.e., inside the past light cone of the point (T, 0)). Numerical evi- 
dence for the existence of a countable family of such solutions was given in [2] 
and the main goal of this paper is to prove this fact rigorously (we note that a 
general variational argument for the existence of infinitely many weak solutions 
of equation ([5]) was given before in [3]). As in [1], we will present two different 
proofs of this result. The first proof, given in Sects. 2 and 3, is rather explicit 
and exploits the conformal invariance of the cubic wave equation is an essential 
way. The second proof, given in the Appendix, is more general in the sense that 
it is based on 'soft' topological arguments. 
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Having established the existence of self-similar solutions in Sect. 3, we analyze 
some of their properties in the second part of the paper. In Sect. 4 we derive some 
remarkable asymptotic scaling formulae for the solutions with many oscillations. 
Sect. 5 is devoted to the linear stability analysis. Finally, in Sect. 6 we show that 
all nonconstant solutions have a singularity outside the past light cone. We point 
out that for the sake of clarity of exposition the paper is written in the 'physics' 
style, however its conversion to the 'epsilon-delta' style is routine and we leave it 
to the mathematically oriented reader. 



2 Dynamical system and local existence 

In the studies of self-similar solutions it is convenient to use hyperbolic polar 
coordinates (s, x) defined by 

T — t = e^^ cosh(x) , r = e~'^ sinh(x) . (6) 

The transformation ([6]) is a conformal transformation of the Minkowski spacetime. 
In the hyperbolic coordinates the Minkowski metric reads 

ds^ = e-^' {-ds^ + dx^ + sinh2(a;) dQ^) , (7) 

where — oo < s < oo, x > and dQ"^ is the round metric on the unit two-sphere. 
The surfaces s = const are hyperboloids with constant scalar curvature — 1 
which foliate the interior of the past light cone of the point (T, 0). Due to the 
conformal symmetry of Eq. the function f{s,x) = ru{t,r) satisfies a simple 
wave equation 

dssf-d,J- . [ =0. (8) 
smh [X) 

The self-similar solutions of Eq. ([2]) correspond to static solutions of Eq. ([8]), i.e., 
solutions f{x) which satisfy the ordinary differential equation (here and in the 
following we denote the derivative by the prime) 

^ sinh^(x) ^ ^ 

on the half-line x > 0. In particular, the constant solution f/o(p) = \/2 of Eq. ([5]) 
corresponds to /o(a;) = v^tanh(x). The remainder of this section and sections 3 
and 4 are devoted to the analysis of solutions of Eq. ([9]). 

In order to obtain a dynamical system formulation we introduce 

Kx)=/(x)-x/'(x), dix) = f'ix). (10) 

Then, Eq. iQ is equivalent to the system of first order equations 

b' = —1 , d' = - , (11) 

sinh^(x) ' sinh^(x) ' 



3 



where f{x) = b{x) + xd{x). Rewriting this system in the form 



/b\' b /6 \^ X* fb \^ 

x[-\ = \ 7T— ( - + d] , xd' = J— i - + d] , (12) 

\x J X sinh (x) \x J sinh (x) / 

and applying Prop. 1 of [1], we infer that there exists a one-parameter family 
{b/x,d){x) of local solutions of Eqs. (fT2l) with boundary condition {b/x,d){0) = 
(0,c), analytic in (x^,c) and defined for all c and |x| < ^(c) with some ^(c) > 0. 
We shall refer to these solutions as 'c-orbits' and denote them by b{c, x) and 
d{c,x). It follows from the above that the c-orbits have the following expansion 
near the origin 

6(c, x) = ^c^x^ + O (x^) , d{c, x) = c - ^cV + O (x^) . (13) 



3 Global existence 

In this section we consider the global behavior of c-orbits. Without loss of gen- 
erality we may assume that c > 0. First, we observe that c-orbits are defined for 
all X > 0, as follows immediately from the existence of the Lyapunov function 

G = 2d^ + — ^ , G' = -2 coth(x)— ^ < . (14) 
smh (xj smh [x) 

Second, c-orbits have simple asymptotic behavior since by f / a/ sinh(x) is 
bounded when x —>■ oo, hence the right hand sides of Eqs. (JTTl) are integrable, 
and therefore b{x) and d{x) have finite limits B and D for x — > oo. Moreover, the 
rapid decrease of the right hand sides of Eqs. ffTTl) implies that the limits B{c) 
and -D(c) are continuous functions of c. 

It follows from the above that f{c,x) ~ B{c) + D{c)x for x oo, hence 
solutions of Eq. corresponding to c-orbits are regular at infinity if and only 
if -D(c) = 0. Now, we will show that there is an infinite sequence c„ with n = 0, 
1, 2, . . . such that -D(c„) = 0. The corresponding globally regular solutions are 
characterized by 6„ = -B(c„) and behave asymptotically as 

fix) =bn- ble-^- + 0{e-'^) . (15) 

Numerically one finds a sequence of such solutions (see Fig.[T]) with n = 0, 1, 2, . . . 
zeros and parameters 6^ ~ c„ ~ (n -t- 1)^ (see Table 1 and Sect. 4). Note that 
for the globally regular solutions the conserved energy functional associated with 
Eq. (ED, 



W)4f (/;./-^)<^. 



(16) 



is finite and, as indicated by numerics, monotonically increasing with n (see 
Table 1). 
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Figure 1: d{c) =D/a vs. b{c) = {B + 8D)/a where a = {{B + SDf + L)2)i/6 fo^. 
positive (sohd) and negative (dashed) values of c. 



n 






E{fn) 


c„ (theory) 


hn (theory) 







V2 


1/3 


1.630626 


1.467029 


1 


9.616283 


-3.578348 


4.62810 


9.991135 


-3.631358 


2 


30.13927 


6.315947 


21.5429 


30.681145 


6.363520 


3 


68.58242 


-9.519976 


64.8053 


69.292246 


-9.563216 


4 


130.5379 


13.13018 


153.071 


131.41603 


13.170001 


5 


221.5967 


-17.10516 


309.116 


222.64408 


-17.142226 


6 


347.3277 


21.41418 


556.682 


348.56798 


21.448919 



Table 1: The parameters of the first few solutions /„ generated numerically and 
their comparison with the asymptotic formulae P3|) and fH4|) . 



Let us introduce the phase function 
/b(x)\ 

<p{x) = arctan -H" , with 0' = > . (17) 

For c-orbits we normalize 0(c, x) by the condition 0(c, 0) = 0. Then 0(c, x) = 
{i — l/2)7r at the i^^ extremum of f{c,x). Furthermore, the limit $(c) = 0(c, oo) 
is a continuous function of c and $(cn) = (n + l/2)7r for a regular solution fn{x) 
with n zeros. 

Integrating Eqs. ffTTl) for c ~ yields 

0(c, x) = + 0{c') , m = c'- + 0{c') . (18) 
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To find the behavior of the phase function for large c we rescale the variables 



F{y) = f{x), y = cx, (19) 

so that Eq. ([9]) becomes 

sinh^(?//c) ^ ^ 

For c ^ oo we get the limiting equation 

F" + ^ = 0, F{y) = y-]:y' + 0{y') , (21) 

whose solution is oscillatory. In terms of the original variables this implies that 
limc_^oo 0(c, x) = oc for any finite x > 0, and hence limc^oo ^(c) = oo. Therefore, 
for each n > there exists at least onci value c„ such that $(c„) = (n + l/2)7r. 
This concludes the proof of existence of a countable family of regular self-similar 
solutions of the cubic wave equation. 



4 Asymptotic formula for 6^ and 

In this section, we use the technique of matched asymptotic expansions to de- 
rive the asymptotic scaling formulae for the parameters of solutions /„(x) in the 
large n limit. The solutions with many zeros are approximately periodic with a 
modulated amplitude. In order to extract the periodic part we factorize / in the 
form 

fix) = aix)vit{x)) , (22) 

with suitable functions a(x) and t(x). Plugging this ansatz into Eq. (jH]) we get 
(denoting t-derivatives by a dot) 

at'^v + (at" + 2a't')v + a"v + — = , (23) 

sinh (x) 

and impose the conditions 

sinh(a;)t' = a , (24) 
2a t' + at" = . (25) 

Differentiating the first equation we get 

a' t" 

_ = _ + coth(a;) , (26) 
a t' 

-'^ Numerics indicates that there for each n there is exactly one c„. 
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which we use to ehminate a from the second one. Thus we obtain 



T = I —41 — = -b(-,-] ^ 4.20654632 . (30) 



^ = -^coth(x), (27) 

and by integration (suppressing an irrelevant integration constant) 

t' = smh-^/^{x) , (28) 

and finally 

Jo sinh^/'^(0 

The length T of the t interval corresponding to < s < oo is 

sinh2/3(^) 2 V6'3 
From Eq. fl24|) we get 

a{x) = sinh^/2(x) . (31) 
Using the expressions for a{x) and t{x) in Eq. fl23l) yields 

v + v^ + hv = , (32) 

with 

3sinh2(x) - 2cosh2(a;) 

= ^ ■ , 2/3/ N • 33 

9smh''/-^(x) 

The behavior of t{x) for a; ^ is t(x) 3x^^^, implying h for x — > 0. 

Introducing t = T — t we find t{x) — * | cosh"^''^(a;) for x — oo and consequently 
h 1/4P. The linear term hv can be neglected except near t = resp. t = 0, 
where it dominates the cubic term. From the boundary conditions for f{x) one 

obtains v{t) ft^ fo^ t ^ resp. v{i) a/I^?^^ for t ^ 0. 

In order to extract the leading behaviour for large b resp. c, we rescale t — *■ 
Yesp. i — ^ i/lP'/^. Furthermore we rescale v hy v ^ c^/^v resp. v iP'^^v. 
In the limit 6, c — oo, neglecting nonleading terms in h we obtain the equations 

w + - = , (34) 

resp. 

i) + + -^t; = , (35) 

for the rescaled variables. The rescaled boundary conditions are v{t) for 
t — i> resp. v{t) \/ 1^^^ fo^ f ^ 0. Numerically one finds that the solutions of 
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Eqs. fl34f35p with these boundary conditions converge very quickly to solutions 
of 

v + v^ = , (36) 

with amplitudes Aq resp. Ai. Numerically one finds A^ 0.90247851 and Ai ^ 
0.82273965. After rescaling v{t) — > v(t/A)/A both solutions tend to the solution 
Fi (t) of Eq. (!36l) with the normalized amplitude and the corresponding period 

r = 4^2 ^^-^ = V2B{\, i) ^ 7.41629871 . (37) 

We fix the phase of Fi such that -Fi(O) = 0. Then the rescaled v's tend to Fi 
with some phaseshifts 9o resp. 9i, i.e. 

v{t/A)/A ^Fi{t + e) for t oo . (38) 

Numerically one finds Oq ^ -1.6225533 and Oi ^ 0.8623512. After shifting the 
two solutions with their respective ^'s they coincide asymptotically. 

Supposing we have a regular solution fn{x) with n ^ 1 zeros, the parameters 
bn and c„ must have been chosen so that the corresponding solutions of Eq. fl32l) 
starting at either end of the interval < x < oo match at some intermediate 
point. Using the discussed asymptotics of the rescaled solutions we obtain two 
conditions matching amplitudes and phases. Matching the amplitudes we get 

c'J'Ao = h^A, . (39) 

To match the phases we assume that we match a solution starting at a; = with 
m zeros with one from x = oo with n — m + 1 at their last zero. The corresponding 
t intervals to and ii must add up to the total t interval T. Taking into account 
the rescalings and the phaseshifts we obtain the condition 

-1/3 /i^^/^ 

Making use of Eq. fl39|l this can be rewritten as 



(n + l)f -(gp + gQ 

Cn ^0 



T , (41) 



which is independent of m and hence from the matching point as required for 
consistency. Thus we obtain the desired asymptotic formula for n ^ 1 



(42) 
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Putting in numbers we get 

3 

(43) 



3.70814935 (n+l) + 0.7602022 



3.7963177 



) 



together with 
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6„ = I ^ I 1.3198462 c„. (44) 



5 Linear stability analysis 

The role of self-similar solutions in the Cauchy evolution depends crucially on 
their stability under small perturbations. To analyze this issue, in this section 
we determine the spectrum of the linearized operator around the self-similar 
solutions fn{x). Substituting the ansatz f{s,x) = fn{x) +w{s,x) into Eq. ([8]) 
and linearizing, we obtain the linear wave equation with a potential 

dssW-d^^W + Vnix)w = 0, Vnix) = - , (45) 

smh [X) 

Since fnix) ~ CnX for x — and /n(oo) = bn, the potential Vn{x) is everywhere 
bounded and decays exponentially at infinity. 

Separating time, w{s,x) = e^^'^C,{x), we get the eigenvalue problem 

L„e = fc^e , L„ = + V;(x) . (46) 

For each n, the operator L„ is self-adjoint on V{Ln) = G i^2[0, oo),^(0) = 0} 
and has a continuous spectrum A;^ > 0. The discrete spectrum depends on n. 
More precisely, we claim that L„ has exactly n + 1 negative eigenvalues. To see 
this, note that Eq. ([8]) is invariant under the following transformation 

s ^ s + -ln (l + 2acosh(a;)e" + a^e^'') , a; ^ tanh"^ ( — | ^ 

2 \cosh(x) + ae*/ 

which is nothing else but the time translation t t + a, expressed in hyperbolic 
coordinates. Hence, each time-independent solution f{x) gives rise to the one- 
parameter family of time-dependent solutions 



/.(.,x) = /ftanh- (^^\] . 
^ ^ V Vcosh(x) + aeV / 



(4J 



The perturbation 5f generated by this symmetry (which we shall refer to as the 
gauge mode) is tangent to the orbit fH8|) at a = 0, that is 

5f = ^^^^^\a=o = sinh(x)/'(x) . (49) 
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Thus, for each n the operator L„ has the eigenvalue fc^ = — 1 associated with the 
eigenfunction ^'■"•'(x) = sinh(x)/4(a;). Since by construction the solution /„(x) 
has n oscillations, it follows that the eigenfunction i'^^\x) has n nodes, which in 
turn implies by the Sturm oscillation theorem that there are exactly n eigenvalues 
below k"^ = —1. This means that the solution fn{x) has at least n unstable modes 
(the gauge mode does not count as a genuine instability). Numerics indicate that 
there are no eigenvalues in the interval —l<k'^< 0, so we claim that the above 
phrase 'at least n' can be replaced by 'exactly n\ however we can prove this claim 
only for the perturbations of the ground state solution fo{x) = v^tanh(x). In 
this case 

Voix) = ^ (50) 

cosh [X) 

is the Poschl- Teller potential for which the whole spectrum is known explicitly. In 
particular, the gauge mode ^^^\x) = sinh(x)/ cosh^(x) is the only eigenfunction. 



6 Behavior of solutions outside the light cone 

The hyperbolic coordinates ([6]) cover only the interior of the past light cone, 
hence in order to see what happens outside the light cone we need to return to 
the similarity variable p and Eq. (Q. The results of Sect. 3 imply that Eq. ([5]) 
has infinitely many solutions which are smooth on the closed interval < p < 1 
and behave as 



and 



[/(p) = c+^(2-cV + (^(p') for p^O, (51) 

f/(p) = 6 + ^(62-2)(p-l) + 0((p-l)2) for p^l. (52) 

We will refer to solutions satisfying the boundary condition ( !52l) as to '6-orbits'. 
Without loss of generality we assume that 6 > 0. Now, we will show that there 
are no 5-orbits which are smooth for all p > 0. The proof consists of two steps. 
In the first step we show that all 6-orbits with < 6 < y/2 become singular at 
p = 0. This will imply that the smooth solutions constructed in Section 3 must 
have bn > for all n > 1. In the second step we show that all 6-orbits with 
b > \pl become singular at some p > 1. 

Stt]) 1. Consider a 6-orbit with < 6 < -\/2 and assume that it exists for all 
p < 1 and is smooth at p = 0. Let us define the function 

Hp)--^- (53) 

Using equation ([5]) we get 

{p-p')U'' + {2-4p')UU'-pU'{2-U') 
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It follows from ([52]) that h{l) = -(2 - b"^) > and h'{l) = -(2 - b'^){b'^ - 4) < 0. 

2 8 
We will show that h'{p) < for all p > and therefore limp^o+ h{p) > h{l) > 0. 

But smoothness at p = requires that h{0) = 0. This contradiction will prove 

Step 1. To show that h'{p) < we show equivalently that the numerator of the 

fraction on the right hand side of Eq. (15^ 

n{p) = (p - p^)U'^ + (2 - 4p'^)UU' - pU^{2 - U^) (55) 

is non-positive. From (l52l) we have n(l) = and n'(l) = -6^(4 — &^)(2 — 6^) > 0, 

hence n(p) is negative in the left neighborhood of p = 1. If n(p) is negative for 
all p < 1 then we are done. Thus, let us suppose that there is a p = i? such that 
n{R) = and n{R) < for R < p < 1. We will show that this is impossible 
because n'{R) > 0. Note that equation 

n{R) = (i? - R^)U''^ + (2 - AR^)UU' - RU'^{2 - f/^) = , (56) 

viewed as a formal quadratic equation for f/', can be satisfied only if the discrim- 
inant (divided by a positive factor 4f/^ for convenience) 

^ = l-R\l-R^){2 + U^) (57) 

is non-negative. A calculation yields 

HI 

n\R) = _ [-(A + R')U' + R'UiU' - 2)] , (58) 

where terms involving U were eliminated using Eq.([56]). If U{R) > V2 then 
the right hand side of Eq. (1581) is manifestly positive. The case < U{R) < V2 
requires more work. In this case Eq. (l56l) has a single negative root 

f 1 - 2i?2 + Ta") u 

Substituting this value into (l58l) we get n'{R) = {positive f actor) ■ N {U, R), where 

N{U, R) = {A + R^){1 - 2R^ + VA) + (1 - R^)R^{U'^ - 2) . (60) 

One can check that in the rectangle < [/ < a/2, < i? < 1 the function N{U, R) 
has no critical points so its minima and maxima occur on the boundary. It is easy 
to verify that N{U, i?) > on all sides of the rectangle, thus N{U, R) > inside 
the rectangle. This implies that n'{R) > and completes the proof of Step 1. 

Step 2: Consider a 6-orbit with b > \/2 and assume that it exists for all p > 1. 
It follows immediately from Eq. ([5]) that f/(p) is monotone increasing. Next, let 

g{p) = p'U{p)\J\p)-^^{\J\p)-2). (61) 
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From ([52]) we get ^(1) = 6(62 -2)/3 ^ ^jg^j^ ^-^at g{p) > for all p > 1. To 

show this let us compute g'{p) at a point where 5f(p) = 0. After a straightforward 
calculation we get 

g'ig = 0) = ((17f/ - 9)p^ + f/ + 3) . (62) 

Since the right hand side of this equation is manifestly positive for p > 1, g{p) 
cannot cross zero from above, hence g{p) > for all p > 1, as claimed. Thus 

U' 1 

>7r^ (63) 



f/2 - 2 6p 

which after integration from 1 to p yields a contradiction 

mjZ^ > ^ exp cm , (64) 

which completes the proof of Step 2. 

It follows from the above analysis that each self-similar solution Un{p) with 
n > 1 is singular at some pn > 1. This fact suggests that these solutions do 
not participate in the Cauchy evolution of smooth initial data which in turn 
corroborates the conjecture that the constant solution Uq = \/2 is a universal 
attractor for blowup solutions (see [6J for what is known rigorously and [7] for 
recent numerical evidence) . 
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Appendix 

We present here an alternative 'soft' topological proof of existence of infinitely 
many smooth self-similar solutions of Eq.(5). 

We already know that given any c there is a unique smooth solution U{p, c) 
of Eq.(5) satisfying f/(0,c) = c defined for all < p < 1. Similarly, there is a 
unique smooth solution f/(p, 6) satisfying U(l, b) = b defined for all < p < 1. 

We express these solutions in terms of polar coordinates, that is we define 



r (p, c) = y/U{p,cy + U'{p,cr , e{p, c) = arctan f 

V U (p, c) 



(65) 
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and similarly 



R{p, h) = ^U{p,hf + U'{p,hY , (3{p, b) = arctan (^^) ■ (66) 

Let po = a/2/3. Since the region {(p, c)|0 < p < PO)C > 0} (resp. {(p, 6)|po < 
p < l,b > 0}) is simply connected, the angle 9{p,c) (resp. (3{p,b)) is defined 
unambiguously once we specify its value at any point in the domain. We set 
^(0, 1) = 0, hence ^(0, c) = for all c> 0. Similarly, we set V2) = 0; then 
-7r/2 < b) < 7r/2 for all b > 0. Next, we define maps 

^:R+3c^ $(c) = {9{po, c), i?(po, c)) e i?^ , (67) 

and 

: i?+ 9 6 ^ ^fe(6) = (/5(po, b) - 2A:7r, i?(po, 6)) G i?^ . (68) 

Note that if ^'fc(^) = $(c) for some b and c, then we have a solution defined on 
the whole interval < p < 1 in the nodal class with index n — 2k (according to 
our terminology from Section 3). 

Lemma 1. limc_>o ^(Po; c) — and limt^o R{Po, b) = 0. 

Proof: Follows immediately from continuous dependence on initial conditions. 

Lemma 2. limc_»oo d{po, c) — —oo and lim6_>oo /3(po, b) — oo. 
Proof: Follows from the asymptotic analysis given in Section 4. 

Lemma 3. For any positive b and c we have 9{p, c) < 7r/2 and /3(p, b) > — 7r/2. 
Proof: We have 9(0, c) = and if 9{p,c) = tt/2, theni U{p,c) = 0,U'{p,c) > 
so 9'{p,c) < 0, contradiction. Similarly, we have /9(1,&) > —n/2 and if P{p,b) — 
-■k/2 then U{p, b) = 0, f/'(p, 6) < so /3'{p, b) > 0, contradiction. 

Lemma 4. If < c < 2, then -7r/2 < 6'(po, c) < 7r/2 and similarly, if < 6 < 2, 
then -7r/2 < /?(po,6) < 7r/2. 
Proof: We define the function 

H{p)^^{l-p')U"-U' + ^U'. (69) 

We have i^'(p) = {3p—2/p)U''^ so i^(p) decreases on (0, po] and increases on [po, 1]. 
If < c < 2 and p < Po then if(p, c) < H{0, c) < 0, hence U{p, c) > (since 
> if ?7 = 0). Similarly, if < 6 < 2 and p > po then H{p, b) < H{1, b) < 0, 
hence C/(p, b) > 0. 

Now we are ready to prove: 

Theorem. For any positive integer n there exist parameters (c^, 6„) such that 
the corresponding solution C/(p, c„) = C/(p, 6„) is in the n*'^ nodal class. 

Proof: If n = 2k then by Lemmas 2 and 3, for any integer A; > 1 we may choose 
cr > cl > 2 (resp. bu > bL > 2) such that 9{po, cl) — —n/2 and 9{po, cr) — 
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-{2k + l)7r (resp. /3(po, h) = 7r/2 and P{po, 6r) = {2k + l)7r). Let cr (resp. 6r) 
be the smallest such c (resp. b). Then — 7r/2 > 6{pq,c) > — (2/c + l)7r for bi < b < 
bji.Then 7r/2 < /3(po, < (2/i; + l)7r for bi <b < bji. Next, we choose m and M 
such that m < r(po, c) < M for cl < c < cj? and m < R{po, b) < M ioi bi < b < 
bR. Finally, by Lemma 1 we choose c < cl (resp. b < bt) such that r(po, c)) = m 
(resp. R{p{), b) — m) and let c (resp. 6) be the largest such c (resp. b). Let Q be 
the rectangle with vertij:^es {-{2k + l)7r, m),(-(2A; + l)7r, M),(7r, m),(7r, Af ). The 
ordered points A = '^{b),B = $(c),C = ^(6r),-D = $(cr) lie on the boundary 
of fl, thus it follows from elementary topology that the curve {$(c)|c < c < Cr} 
from B to D and the curve {"^ k{b)\b < b < bR} from >1 to C must intersect. 

If n = 2A;+ 1 we can repeat the above argument with 6 < making the obvious 
modifications in Lemma 3 (/3(p, b) > 7r/2) and Lemma 4 (7r/2 < /5(p, b) < 37r/2). 

References 

[1] p. Bizoh, D. Maison, and A. Wasserman, Nonlinearity 20, 2061 (2007) 

[2] P. Bizoh, T. Chmaj, and Z. Tabor, Nonlinearity 17, 2187 (2004) 

[3] V.A. Galaktionov and S.I. Pohozaev, Quart. Appl. Math. 61, 583 (2003) 

[4] P. Breitenlohner, P. Forgacs, and D. Maison, Commun. Math. Phys. 163, 

141 (1994) 

[5] O. Kavian and F. B. Weissler, Commun. PDE 15, 1421 (1990) 
[6] F. Merle and H. Zaag, Math. Annalen 331, 395 (2005) 
[7] P. Bizoh and A. Zenginoglu, math- ph/081 1.3966 



14 



